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Abstract 

Fluctuations associated with stretched E2 transitions from high spin levels in nuclei 
around ^^^Yb are investigated by a cranked shell model extended to include residual 
two-body interactions. It is found that the gamma-ray energies behave like random 
variables and the energy spectra show the Poisson fluctuation, in the cranked mean field 
model without the residual interaction. With two-body residual interaction included, 
discrete transition pattern with unmixed rotational bands is still valid up to around 
600 keV above yrast, in good agreement with experiments. At higher excitation energy, 
a gradual onset of rotational damping emerges. At 1.8 MeV above yrast, complete 
damping is observed with GOE type fluctuations for both energy levels and transition 
strengths (Porter-Thomas fluctuations). 



1. Introduction 

Recently, improved detectors and analysis techniques have focused the attention 
on the off-yrast nuclear structure of deformed nuclei at high spin. For example, it is 
expected that the rotational band structure typical of the near-yrast regions is smeared 
out at an intrinsic excitation energy of order of £^a; ^ 1 MeV measured from the yrast 
in rare-earth nuclei There are also arguments that this damping property of the 
rotational motion is related to onset of chaotic behavior in deformed nuclei ^'^'^) . 

Information of the nuclear structure of rotating nuclei is usually obtained from 
identifying the E2 transitions in rotational bands near yrast, while the gamma-rays 
emitted from levels at higher excitation energy form the so-called "quasi-continuum" 
spectra in which individual transitions are not distinguished because of the high level 
density and finite detector resolution. The fluctuation analysis technique, which has 
been developed recently is useful for a study of the quasi-continuum spectra. From 
the analysis of fluctuations in the ridge region of 2-D E.y-Ej spectra, an effective number 
of the E2-decay paths associated with decay along rotational bands is extracted. The 
fluctuation analysis technique is based on the assumption that the gamma-ray energies of 
the rotational E2 transitions behave like random variables, as if they are like rain drops 
falling at arbitrary positions ^). This is the key assumption in making a connection 
between the fluctuation of counts in the quasi-continuum spectra and the number of 
rotational bands passed in the decay. 

It should be emphasized that we actually do not know the underlying fluctuation 
behavior of the E2 transitions. Consider rotational bands in the near-yrast region which 
do not interact strongly with each other. One may expect that if some dynamical 
quantum numbers or symmetries persist, they imply regular patterns or correlations in 
the rotational transition energies. On the other hand, the rotational perturbation at 
high spin may be able to destroy such regularity by breaking the dynamical symmetries. 
Hence the validity of assumption of random E2 transition energies, adopted by the 
fluctuation analysis, may depend on whether any dynamical symmetries in the rotating 
nuclei survive. 

The E2 transitions are expected to show very different fluctuation behavior when 
the rotational damping sets in. In this case, the E2 transitions from an initial level 
fragment into many final levels because many rotational bands with different intrinsic 
configurations are strongly mixed. One may readily expect that the wave functions of 
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the strongly mixed levels are very complicated, and that strengths of individual com- 
ponents of the fragmented E2 transitions strongly fluctuate, reflecting the complexity 
of the mixed wave functions. A description of the strength fluctuation in the strongly 
mixed levels may be governed by the random matrix model (the Gaussian Orthogonal 
Ensemble), in which the strength fluctuations show the Porter-Thomas distribution ^). 
It is not clear, however, whether the Porter-Thomas strength fluctuation is realized in 
the E2 transition matrix elements in the region in interest, with intrinsic excitation 
energy of a few MeV. 

In this paper we examine theoretically the fluctuations associated with the E2 tran- 
sitions in rotating nuclei on the basis of a model for the high spin off-yrast levels. In 
particular we shall investigate whether randomness of the transition energies is realized 
for the non-interacting rotational bands, and whether the strengths in the damped E2 
transitions exhibit Porter-Thomas fluctuations. 

2. Extended cranked shell model 

The model we use is an extension of the cranked shell model, in which we take 
into account not only the cranked single-particle Hamiltonian but also the two-body 
residual interactions which cause the rotational damping. As the single-particle fields 
we adopt the Nilsson potential and the cranking term —u)Jx, but neglect the static 
pairing gaps, since the levels emitting the quasi-continuum gamma-rays are in the high 
angular momentum region around 30^ to 60^, where the pairing gap is expected to be 
negligible. 

It is straightforward to introduce the independent-particle configurations using 
single-particle orbits in the cranked Nilsson field. Many-particle and many-hole (np-nh) 
excitations should be included to form a suitable basis for the off-yrast levels. The 
independent-particle configuration is hereafter denoted as |/i(a;)). As an ansatz, it is 
assumed that the np-nh configurations correspond to the rotational bands which would 
exists in the off-yrast region if there were no interactions between the rotational bands. 
A residual two-body interaction Hres is assumed between the rotational bands. Thus 
the Hamiltonian in this model is expressed: 

H = hNiisson — '^Jx + Hres (1) 

as a function of the rotational frequency uj (using the convention h — 1). Details of the 
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residual interaction are discussed in sect. 4. 

The Hamiltonian H is diagonalized within a suitably truncated basis set of {|yu(u;))}. 
The diagonalization is separately made among the basis states with the same parity tt 
and signature a at fixed uj. The resultant levels are denoted by i, their energy by 
Ei{ijj) and their wave functions by \i{oj)). The routhians of the levels are defined by 
e\{uj) — Ei{uj) — Ej.f.f{ijj), where the reference energy Ej.f,f{ijj) is chosen as the lowest 
energy level in the numerical calculations. 

The rotational stretched E2 transitions are calculated as follows. As far as we con- 
sider well-deformed nuclei with stable shape, it may be assumed as a good approximation 
that intra-band transitions dominate if no interaction between the bands exists and the 
quadrupole moments are the same (constant Qq) for all the bands. The E2 transition 
strength and the associated 7-ray energy can be calculated as 

Sij = gQl{3{uj-2/J)\t{uj)f , (2) 

E^^2u; + e'iiu) - e'j{cu - 2 /J) = Eij , (3) 

for a transition between levels i and j. Here \j{u) — 2/ J)) is a. solution of the extended 
cranked shell model at the rotational frequency uj—2j J ^ J denoting the average moment 
of inertia. In this way, i and j represent states at angular momentum / and 1 — 2, 
respectively. The factor g in Eq.(2) represents the geometrical factor which arises from 
the rotational I?- function. Only the overlap matrix element — 2/J)\i{ijj)) plays a 
role in the following analysis since we always normalize the E2 strengths. 

The above framework is similar to the one employed by S. Aberg '^''*) in many re- 
spects. References ^'^), however, introduce an additional procedure to transform the 
energy of the np-nh basis \iJi{oj)) to the lab-frame energy as a function of the angular 
momentum /, and the diagonalization is made at a fixed /. We do not follow this proce- 
dure here because it is not expected to influence the statistical properties of the results, 
although individual levels may differ in both descriptions. The procedure in refs.^'^) may 
be useful for direct comparison with experiments, especially when describing E^ — Ej 
spectra. The procedure described above is a more straightforward theoretical framework 
and simpler in practical numerical calculations. 
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With the residual interactions incorporated, the extended cranked shell model can 
describe both the rotational bands near the yrast line and the levels far off yrast, where 
strong bands mixing may occur. Nevertheless, we shall divide our analysis into two parts 
for convenience of discussion. We first neglect the residual interaction, as in the con- 
ventional cranked shell model. The intra-band E2 transitions in the ncar-yrast regions 
are then described by this approximation since it is expected that the rotational bands 
in this region do not mix with each other. We shall discuss the statistical properties for 
those intra-band rotational transitions in sect. 3. In sect. 4, we discuss the results of 
the extended cranked shell model where the residual interactions are incorporated, con- 
centrating on the E2 transitions in the far-off-yrast region where the rotational damping 
is found. In sect. 5 a brief discussion of the transitional region between unmixed and 
mixed rotational bands is given. 

3. Non-interacting rotational bands 

In the cranked shell model without residual interaction, the rotational bands are 
simply described as independent-particle configurations assigned to the cranked single- 
particle orbits (called np-nh bands hereafter) and only the intra-band stretched E2 
transition is considered. The transition strength is almost constant, irrespective of the 
configurations of the bands, provided that the deformation is stable. On the other hand, 
the transition energy varies from band to band because of differences in the alignment. 
Thus it is appropriate to focus our attention on the fluctuation property of the transition 
energies. 

3.1. DISTRIBUTION OF 7-RAY ENERGIES 

The gamma-ray energy of the the stretched intra-band E2 transitions of a np-nh 
band is expressed as 

= 2a;^(/) 

2 (4) 
^2a; + -(7-J») 

in terms of the expectation value J^{oj) of Jx for an independent-particle configuration 
at a;. In expressing the rotational frequency oJ^il) which corresponds to spin value 
7, we make the approximation as 7 = J fj_{uj ~ Jn{i^) + J{oJfj,{I) — oj) with J 
being the moment of inertia. Eq.(4) can also be derived from Eq.(3). It should be 
noticed that the stretched E2 transitions occur successively following the spin sequence 
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•••7 + 2— — 2---. In order to investigate the properties of spectra containing 
all these sequential transitions we displace the spectrum by ±4/ J", ±8/^7, etc, and 
fold all the displaced spectra. Here the energy unit 4/ J" of the displacement is assumed 
to be 0.06 MeV which is appropriate for well-deformcd nuclei with A ~ 170. The folded 
Ej spectra are equivalent to the spectra of J^{ui) modulo 2, which we denote by A J^. A 
typical channel width of the experimental spectra is about one tenth of the 7-ray energy 
unit 4/ J'. Correspondingly, we calculate the distribution of AJ^ = modulo 2 

after sorting the A J^'s into 10 channels. Figure 1 shows a typical result for ^^^Yb. The 
Nilsson parameters in ref.^) and the deformation parameters €2 — 0.255, €4 — 0.014 in 
ref.^) are used. Here the spectra in fig.l are calculated for the lowest 20 bands (fig. 1(a)) 
with the signature and parity quantum numbers (+, 0) as well as for the next 20 bands, 
the third 20 bands and so on, altogether 10 bins of bands. 

Fig.l 



3.2. FLUCTUATIONS OF DISCRETE 7-RAY SPECTRA 

The spectra in fig.l display a quite random fluctuation. In order to analyze the 
fluctuation quantitatively, the second moment or the variance of the spectra is calculated 
in accord with the fluctuation analysis of the experimental data ^). The second moment 
//2 here is defined by fi2 = CEn'^l)/ ^channel - {{Y^n^n) /Nchannelf where rUn is the 
number of AJy^'s in the n-th channel. Prom the obtained for the 10 bins we calculate 
the average value at different rotational frequencies u. The result is plotted in fig. 2 as 
a function of uo. For comparison we show in the same figure the expected limit if 
the distribution of AJ^ is completely random, i.e. obeys a Poisson distribution. The 
random limit is calculated by assuming that the distribution of number rrin in channel 
n is binomial. The resulting average value of //2 is (1 — ^ channel) i^'^ = 1.8 with the 
mean number //i = 2, and the statistical sampling allows a deviation within 0.26. The 
standard deviation of 1x2 in the limit is cr(/i.2) = \/ '2, /N channel A*i — 0.89. The random 
limit is indicated in fig. 2 with the horizontal line and the expected deviation with the 
dashed lines. 

Pig.2 
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It is clearly seen in fig. 2 that the fiuctuation in the calculated spectra follows the 
random limit for rotational frequencies oj higher than 0.3 MeV. It is also seen that there 
is a systematic deviation from the limit for low rotational frequencies. We find that this 
deviation arises from the fact that time-reversal doublets in the cranked single-particle 
orbits are less dissolved for lower cu, because the rotational perturbation is weak. This 
causes a degeneracy in AJ^, which results in larger fluctuations. It should be noted, 
however, that the cranked shell model does not take into account effects which influence 
the spectra at low spin. For example, the intrinsic angular momentum components 
along the symmetry axis {K components) is expected to split the degenerate states 
which have different K values. The effect of the K component on the 7-ray energy is 
evaluated according to Eq.(A.5) given in the appendix. We estimate values of K for 
rotational bands by assigning them to those that are obtained at a; = 0. This gives a 
order-of- magnitude estimate of the K effect. The resulting second moments of the 7-ray 
energy spectra arc shown in flg.2 as triangles. One can see that the inclusion of the K 
effect enhances the tendency toward the random limit. We conclude that the random 
limit is achieved for the high spin region a; ^ 0.3 MeV or ^0.6 MeV. 

It is neccesary to check the above conclusion because our sampling of the rotational 
bands in the above analysis may be somewhat artiflcial. We sample the lowest 200 np- 
nh bands for each (tt, a) while the number of the unmixed rotational bands is order of 
a few tens according to the experimental fluctuation analysis ^) . We therefore perform 
a physically more relevant analysis by sampling only the bands near yrast. Prom np-nh 
bands with all combinations of (7r,Q;), the lowest 20 bands are picked up. Associated 
AJ^/s are sorted into the channels in the same manner as in flg.l. To obtain a set 
of the second moments /i2 of spectra, we perform the same calculation for 10 even- 
even nuclei around ^^®Yb, namely, 162-I66gj,^l64-170Y13^ and ^^^~^'^^Hf. From the set 
of yU2's, we calculate the mean value and the standard deviation (7(/X2) as a function 
of uj. The result is shown in flg.3, also with the K effect taken into account. From 
the comparison with the average //2 and (7(/i2) of the random Poisson distribution, the 
following picture emerges: The. intra-band ^y-ray energies of the rotational bands in the 
independent-particle cranked shell model appear to be randomly distributed, at least at 
high spin. 
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Fig.3 



4. Mixed rotational bands 

4.1. RESIDUAL INTERACTION 

The two-body residual interaction consists of the monopole and quadrupole pairing 
interaction and the quadrupole-quadrupole interaction, which represent large compo- 
nents of the two-body matrix elements. We also introduce additional components of 
the two-body matrix elements, which may arise from the other multipolarities. By 
assuming the additional components are structureless, we give them simply in terms 
of random numbers following the Gaussian distribution around zero mean value. The 
strength of the residual interaction is fixed as follows. For the monopole pairing in- 
teraction we adopt a standard value Gq = 20 /A MeV. We fix the quadrupole pairing 
strength by G2/G0 = 35/3i?^ which is derived from the multipole expansion of the 
delta-force^). The quadrupole-quadrupole interaction strength is given by the selfcon- 
sistent values for both iso-scalar and iso-vector components with the polarization 
effect {xnn = Xpp = 0-7Xself, Xnp = ^-^Xself where the Xself is the selfconsistent value). 
In the numerical calculation, we used the stretched quadrupole operator (instead of the 
ordinary ones) for practical reasons. 

The strength of the random two-body matrix elements is fixed as in ref.^), namely the 
spreading width = 0.039 (y^) E^^"^ MeV of independent-particle configuration is 
extracted from the systematics of the single-particle spreading width, and related to the 
random two-body interaction through the Fermi golden rule vjj, — 27rp2-body'^^^- The 
density P2-body of levels interacting through the two-body interactions is calculated 
explicitly in the cranking model without the interactions and it is parametrized as 
P2-body — 22£'^/^ MeV~^. The root- mean-square value v of the random two-body 
matrix elements is thus fixed to f = 16 keV. 

It can be argued that there is an uncertainty in the residual interaction strength. 
However this uncertainty is not crucial in the following analysis as long as we do not 
aim at very accurate numbers in the excitation energy. As is discussed in ref. ^), the 
two-body interaction strength scales proportionally to the spacing between interacting 
configurations, i.e. (i2-body = l/P2-body oc E~^/'^. If the scaling is a good approxima- 
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tion, a factor 2, for example, in the interaction strength results in a scahng of 2^/3 .^1.6 
in the excitation energy. This should, however, only be taken as a guideline since we 
are investigating properties concerning the full set of states. 

4.2. CHOICE OF BASIS 

In diagonalizing the Hamiltonian matrix the basis of np-nh configurations is trun- 
cated in the following way. First, we construct all the np-nh configurations that lie 
below Ex = 3.0 MeV above the lowest energy configuration. Next the diagonal energies 
due to the residual interaction arc calculated for these configurations. The basis of the 
diagonalization consists of the lowest configurations including the diagonal interaction 
up to a certain number, eg, 600, 1000, or 1400, for each (tt, a). The levels at a; — 2/ J" 
are calculated using the same basis as defined at u;. Figure 4 shows the density of levels 
calculated thus for the levels with (— , 1) in -"^^^Yb at a; = 0.5 MeV. Results with the 
three different truncations are plotted to display the convergence. The basis of 1000 
states is found to be sufficient for describing the levels up to 2.0 MeV from the yrast, 
amounting to over 200 levels for each (n.a). In the following the calculation with 1000 
states is discussed. Most of the figures presented below are for calculation on ^^^Yb 
at u! = 0.5 MeV, although a few include the results from the 10 neighboring even-even 

nuclei 162-166Er,164-170Yb, and 168-172Hf. 

Fig.4 



4.3. TRANSITION STRENGTH FLUCTUATIONS 

When the residual interaction and band mixing are included in the calculations, we 
see that rotational transitions from each level branch out into many final levels. A clear 
evidence of this manifestation of rotational damping is shown in fig. 5, which depicts 
four typical examples of the calculated rotational strength distributions plotted as a 
function of 7-ray energy. The strengths shown in each quarter correspond to decay from 
individual levels. The four levels with (tt, a) — (— , 1) from different regions of excitation 
energy demonstrate how the rotational decay will change its structure as the excitation 
energy is increased. 
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Fig.5 



For the yrast level, the rotational strength concentrates into almost only one tran- 
sition, indicating the dominance of the rotational band structure. For the 6-th level 
at Ex = 758 keV, the strength fragments into several pieces. We may consider this 
fragmentation as damping of the rotational transitions, although this is rather subtle 
because of the small number of fragments. The situation becomes clearer for the 24-th 
level {Ex = 1.133 MeV) and for the 200-th level {E^ = 1.948 MeV): the rotational 
strength from these levels splits into a large number of fragments and there appears a 
smooth profile from which a damping width can be extracted. Thus we see in fig. 5 the 
onset of rotational damping around several hundred keV. 

Next, let us focus our attention on the statistical properties of the rotational tran- 
sitions, in particular, on those associated with the rotational damping. As seen in fig.5, 
the height of the peaks or magnitude of the rotational strengths fluctuates significantly 
around the smooth profile of the strength distributions. This fluctuation of the strength 
is what would be expected for levels with complex mixing. To clarify the nature of the 
strength fluctuation, we shall discuss it on a quantitative basis. In particular, we shall 
discuss whether it exhibits the Porter- Thomas fluctuations or not. 

In analyzing the fluctuation in the rotational strength it should be noted that the 
strength distribution exhibits an overall smooth proflle, which may be characterized by 
the central peak position and the rotational damping width (or FWHM) . The strength 
fluctuation seen in flg.5, therefore, should be measured relative to the smooth proflle. 
It is appropriate to single out the relative fluctuation by normalizing the E2 strengths 
with respect to the smooth proflle. For a transition from itoj, the normalized strength 
is deflned by 

^ ' 

where (Sij) represents the strength in the absence of fluctuations. By definition, the 
normalized strength should have unit mean. The smoothed strength {Sij) should be 
a function of the 7-ray energy Eij and the excitation energy of the levels. One can 
evaluate it by 

{S.,) = S{E^j)/p{e'j) , (6) 
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where S{E) is the smooth profile of the rotational strength distribution function. The 
smooth strength function S{E) is extracted by taking a Strutinsky-type smoothing of 
the microscopically calculated strength distribution 

ij i 

where the average over initial levels yields a reliable smoothing. In the Strutinsky 
smoothing, a Gaussian width of 50 keV is used. The level density p(e) is determined by 
a fitting to the constant temperature formula^^) p(e) = T~^e^^~^°^/'^ . 

Fig.6 



In order to sec the excitation energy dependence, the initial levels with the same 
(tt, a) are grouped into bins, each of which contains 50 levels. The excitation energies 
of these bins are approximately 0.0 — 1.4 MeV for 1 — 50-th levels, 1.4 — 1.7 MeV 
for 51 - 100-th, 1.8 MeV for 101 - 150-th, and 1.9 MeV for 151 - 200-th. The 
rotational strength distribution calculated for these bins are shown in fig.6 as well as 
the smoothed strength function S{E) plotted with the solid curves. We sample the 
normalized strengths of the transitions whose 7-ray energy Eij lies within the interval 
of the full width at the one eighth of maximum of S{E). For the lowest bin, for example, 
those transitions picked up satisfy 881 keV < E^j < 1170 keV. 

Fig.7 



After the sampling we calculate distribution of the normalized strength Sij. The 
result is shown in fig. 7 , together with a Porter-Thomas distribution 

P{s)ds = (27rs)-V2e-V2rfs (8) 

for comparison. It is clear from fig.7 that that the distribution of the normalized 
strengths approaches that of Porter- Thomas for increasing excitation energy and al- 
most reaches it in the region E^ ^ l.S MeV. Deviation from the Porter- Thomas is seen 
in the lowest bin which corresponds to Ex — 0.0 — 1.4 MeV (also slightly in the next 
bin). 
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Fig.8 



In the above analysis it has imphcitly been assumed that the property of the strength 
fluctuation does not depend on the 7-ray energy. Let us show this for the highest bin 
which clearly displays Porter-Thomas fluctuations. For this purpose we subdivide the 
samples of the normalized strength into two groups according to 7-ray energy. One 
subgroup is for the transitions whose 7-ray energies are close to the central peak of 
the rotational strength function. In practice we choose the transition energies within 
the Full- Width-Half-Maximum of the smoothed strength function. The other is for 
transition energies which lie in the tail portion of the strength function (between the 
half and one eighth of the maximum value). Figure 8 shows the result of this analysis 
for transitions in the fourth bin (151 - 200-th levels at ~ 1.9 MeV in ^^^Yb) of 
figs. 6 and 7. Regardless of whether the transitions arc from the central or from the tail 
portions, the strength fluctuation is found to obey the Porter-Thomas distribution. 

It should be noticed here that the Porter-Thomas distribution manifest itself only 
for the relative strengths which are normalized with respect to the smooth profile of the 
rotational strength function. This explains why the Porter-Thomas distribution emerges 
here while it does not in the analysis made in ref.^). 

4.4. BRANCHING NUMBER 

Let us analyze the strength fluctuation again, but this time by looking into it in terms 
of a quantity which is closely related to the experimental observables. In the experimen- 
tal fluctuation analysis ^) the number of the rotational decay paths n^p^^j^ is extracted 
from the second moment of the spectral fluctuation. Assuming one-dimensional spectra 
(or observing one-step transitions), the number of paths is given by 



where fi represents the probability for flow passing through the level i, and Wij the 
probability for the rotational transition to pass from level i to level j (with normalization 




(9) 




(10) 
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^jWij — 1). The quantity nbranch(0 ^ measure of the number of transitions which 
branch out of level i (or simply the branching number). With the microscopic model 
the branching numbers are determined by the transition probability as given by 



(11) 



j 



where the energy dependent factor oc is neglected. On the other hand, the average 
branching number expected for a group of levels whose E2 strength displays Porter- 
Thomas fluctuations can be estimated analytically by using a suitable profile of the 
rotational strength function. In fact, the Gaussian form with FWHM of the rotational 
damping width Trot is adopted. The average branching number in the Porter- Thomas 
limit is then given by 



Figure 9 compares the Porter- Thomas hmit Eq.(12) with the average branching 
number calculated microscopically by means of Eq.(lO). In Eq.(12) we adopt Trot — 200 
keV. This value agrees well with the FWHM's extracted from the smoothed strength 
function shown in fig.6. Actually they are approximately 170 keV and the agreement is 
fulfilled within an accuracy of 15 %. The vahdity of Eq. (12) supports again that the 
Porter-Thomas fluctuation is realized in the damped rotational transitions. 

4.5. ENERGY LEVEL FLUCTUATIONS 

It may be interesting to look into the fluctuation in the energy levels because it 
is expected to show the properties of the GOE random matrix theory if the Porter- 
Thomas strength fluctuation is a consequence of the complex mixing of the rotational 
bands. Here we test the A3 statistics ^) as a measure of the energy level fluctuations. 
The calculated spectra are unfolded for each (tt, a) by use of a fltted level density, 
parametrized according to the constant temperature formula. An example of the fltting 
is shown in the inset of flg. 10. The mean value of A3 is calculated as a function of 




(12) 



0.5 Trot p{Ei) 



Fig.9 
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the interval length L (the unit is the average spacing). As in the preceding section, we 
subdivided the levels into bins of 50 levels. 

Fig.lO 



Figure 10 shows that the A3 statistics approaches the GOE hmit with increasing 
excitation energy and that the limit is reached for the levels above #150 at £'3; ~ 1.8 
MeV. It is remarkable, compared to the results shown in fig. 7, that the tendency toward 
the GOE limit is completely consistent with that of the strength fluctuation. Thus it 

can be concluded, from both the strength fluctuations and the energy level fluctuations, 
that the mixing of the rotational bands is as complex and random in the energy region 
Ex ^ 1.8 MeV as described by the GOE random matrix theory. Concerning the energy 
level fluctuations, the results are consistent with the work by S. Aberg ^'^). 

We would like to emphasize the important implication of the above result for the 
rotational transitions in the damping region. A spectrum which shows GOE fluctuations 
is known to be "rigid" with very small fluctuations. As a consequence, the fluctuations 

in transition energy play a relatively unimportant role in the damping region, while the 
strength fluctuation becomes the dominant origin of the spectral fluctuations ^^). 

5. Onset of damping 

As discussed in connection with flg.5, the model predicts both the damped tran- 
sitions and the discrete rotational transitions which characterize unmixed rotational 
bands near yrast. Thus it allows us to investigate the onset of the rotational damping 
and also to analyze fluctuation properties of the discrete rotational transitions, which 
was discussed in the preceding section without including the residual interaction. 

Fig.ll 



5.1. UNMIXED ROTATIONAL BANDS 

Figure 11 illustrates the behavior of rotational transitions in the energy region where 
the damping is about to set in. The branching numbers calculated for individual levels 
are also shown in the flgure. A quantitative criterion for deflning levels associated with 
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unmixed rotational bands may be given by a condition that the branching number of a 
level is less than 2. Accordingly the onset of the damping will be marked by ribranch > 2. 
If this criterion is assumed we count 22 unmixed rotational bands in ^^^Yb, and 30.9 
bands on average for the 10 rare-earth nuclei around spin 42^. Considering that the 
pure rotational bands carry essentially 100 % strength of the intra-band transition, one 
may also define unmixed rotational-band levels as those whose strongest transition has 
a large fraction of the total strength, e.g. l/\/2 = 71%. This gives us gives us 19 
rotational bands for ^^^Yb and 28.3 if averaged over the 10 neighboring nuclei. The two 
criteria give almost the same number. The above number of the unmixed rotational- 
band levels is consistent with the experimental fluctuation analysis which reports 
approximately 30 bands in ^^^Yb. 

Since the residual interaction effect was not included in the preceding analysis using 
the independent-particle cranked shell model in sect. 3, we re-examine the fluctuation 
property of the transition energies associated with the discrete rotational transitions 
defined by the above criterion. For this purpose, we select the lowest 20 levels which 
satisfy ribranch < 2 and use their strongest transitions. The second moment 112 which 
represents the fiuctuation of the selected transition energies is calculated for each nucleus 
at (jj = 0.5 MeV in the same way as in sect. 3. The average value of /X2's over the 10 
nuclei is found to be 2.16, and their standard deviation is 0.80. These numbers compare 
well with the random Poisson limits, 1.80 and 0.89, respectively. Thus the residual 
interaction does not modify the conclusion of sect. 3, concerning the randomness of 
transition energies. 

5.2. TRANSITION STRENGTH FLUCTUATIONS 

The energy for the onset of the rotational damping may be defined although the 
damping sets in rather gradually with increasing excitation energy as seen in fig. 11. 
Prom fig. 9, which shows average excitation energy dependence of the branching number, 
and with the criterion ^branch > 2, the onset energy for l^Syi^ be estimated to be 
approximately Ex ~ 600 keV above yrast. 

Above the onset energy, the transitions from each level begin to branch into many 
small pieces and strong fiuctuations in the strengths of the fragmented transitions set 
in. Although in fig.7 almost pure Porter-Thomas fiuctuations are found in the energy 
region Ex ^ 1.8 MeV, a clear deviation from Porter-Thomas is found for lower exci- 
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tation energies. To describe the energy dependence of the strength fluctuation more 
quantitatively, the lowest 50 levels are subdivided into smaller bins of 10 levels and the 
distribution of the normalized strengths is calculated for each bin, producing the result 
shown in fig. 12. In this figure, results for 10 nuclei are included to obtain good enough 
statistics. It is seen in fig. 12 that deviations from Porter-Thomas fiuctuations are more 
significant at lower excitation energies. Strong deviations are found not only in the 
lowest bin {Ex — 0.0 — 0.9 MeV) which contains the unmixed rotational-band levels but 
also in the other bins where the rotational damping becomes dominating. This indicates 
the presence of a transient region between 0.6 MeV and 1.8 MeV where the rotational 
transitions are already fragmented significantly, while the strength fluctuations have not 
yet reached the Porter-Thomas limit. 

Fig.l2 



5.3. LEVEL SPACING STATISTICS FOR UNMIXED BANDS 

From the viewpoint of discrete spectroscopy ^^), fluctuation properties in the energy 
levels near the yrast line are of interest, although the spin region in the present analysis 
is higher than current discrete spectroscopies can reach. Figure 13 displays the nearest 
neighbor level spacing distribution for the rotational-band levels near yrast. As a ref- 
erence, a result for cranked shell model without the residual interaction is also shown. 
In this flgure only the lowest energy levels satisfying ribranch < 2 (the lowest 5 levels for 
the independent-particle model) are sampled for each (tt, a) from the 10 nuclei. This 
gives 121 spacings in the case of two-body interaction included, and 160 spacings for 
the independent-particle model. The energy dependence of the fltted level density is 
taken into account in normalizing the spacings (cf. sect. 4. 5 and fig. 10). The result for 
the independent particle model resembles the Poisson distribution. With the residual 
interaction included, on the other hand, there are less small spacings than expected 
from the Poisson distribution. The level spacing distribution appears to lie in between 
the Poisson and the GOE hmits. 

Fig. 13 
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6. Conclusion 

The E2 transitions in the off-yrast region of well-deformed rare-earth nuclei are 
investigated on the basis of an extended cranked shell model which describes the np- 
nh bands as well as the mixing among them. Using the cranked shell model without 
the residual interaction, the fluctuation in the gamma-ray energies of the rotational 
transitions associated with the near-yrast non-interacting rotational bands are discussed. 
It is found that the gamma-ray energy spectra at high spin show Poisson fluctuations 
as if the gamma-ray energies were random variables. The model with the residual 
interaction is found to give rise to the damping of the rotational transition above an 
excitation energy {Ex ~ 600 keV). It is seen that the magnitude of the E2 strength 
fluctuates for the fragmented components of the damped rotational transitions. Further 
up in excitation energy {Ex ^1.8 McV), the strength fluctuations arc found to obey the 
Porter-Thomas distribution. In the same energy region, it is also found that the energy 
level fluctuations distribute according to the GOE random matrix model implying that 
the mixing of rotational bands is as complex as in the random matrix theory. Using 
the properties of the random matrix theory as a criterion for the chaotic behavior of 
quantum systems ^^), this region appears to be chaotic. The energy region just above 
the onset of damping appears as a transient region where the strength fluctuation show 
significant deviations from the Porter-Thomas distribution. 

In the model we adopted, the boundary energy for the onset of rotational damping 
is about 600 keV, and the boundary for sharing full Porter- Thomas strength fluctuation 
is around 1.8 MeV. These numbers, however, should be considered with reservation 
since they may depend on the residual interaction as well as its strength. It is therefore 
important to extract these energies in more detailed experiments. 

On the theoretical side, it seems interesting to study how the situation changes 
when more realistic residual interactions are used. Another interesting aspect of band 
mixing concerns the correlations in the wave functions, which extend over several steps 
in angular momentum. These correlations are relevant for a more precise comparison to 
two and higher fold experimental spectra. Such investigations are in progress. 
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Appendix 

With the presence of a K component of the angular momentum, the rotational 
frequency along the rotational axis can be defined by 

, , E(I + 1) - E(I ^ 1) , , , 



IM) = -JP - ■ (A2) 

For a np-nh band in the extended cranked shell model, in which the angular momentum 
Ix along the rotational axis is calculated as a function of uj &s Ix = Ji^{i^), Eq.(A.2) is 
used to determine cu which corresponds to the angular momentum I. 

^JI^ - Kl = J^K(/)) ^ J^(a;) + J{u^{X) - u), (A.3) 

where is the K value for the band and a; is a reference frequency. Eq.(A.l) reduces 
to an equation which gives the associated 7-ray energy, 

81 I 
E,^2uj^{I)^ = 2u,{I) = . (A4) 

V - 

Combining (A. 3) and (A. 4), we get 




(^.5) 



In deriving the final expression, we neglect the /-dependence in the geometrical factor 
// \J ~ ^l by replacing it with Jui^Y + ^fil ^tJ^i^)- Eq-(A.5) can be used in place 
of Eq.(4). 
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Figure captions 

Fig.l The spectra of AJ^ = modulo 2 of the basis np-nh configurations for (tt, a) = 
(+, 0) bands in ^^®Yb at a; = 0.5 MeV. The first bin contains the first 20 bands in 
order of the routhian energies, the second bin is for the next 20 bands, and so on. 

Fig. 2 The averaged second moment //2 of fluctuation defined for the A spectra of the 
np-nh basis bands with (+,0) in ^^^Yb as a function of the rotational frequency cu, 
plotted with closed circles. The value at a; = 0.5 MeV corresponds to the spectra 
shown in Fig.l. The triangles arc calculated for the spectra in which the effect 
of the K components on the 7-ray energy is taken into account. The random 
(Poisson) limit is indicated by the horizontal line. The dashed lines represents the 
uncertainty expected for the average 112 in the random limit. 

Fig. 3 The second moment 112 of fluctuation in AJ^ spectra for the near-yrast np-nh 
bands. The lowest twenty bands are sampled from each nucleus. The average 
value of the /X2's calculated for the ten nuclei (sec text) is plotted with closed 
circles. The bars represent the standard deviation of fi2^s. The random limit for 
the average ^2 is indicated by the solid line and the thick bar indicate the limit 
standard deviation. 

Fig.4 The density of levels calculated in the extended cranked shell model for the levels 
with (— , 1) in ^^^Yb at a; = 0.5 MeV. Results with three different truncations 

(number of basis states 600, 1000, and 1400) are plotted. The Fermi-gas formula 
with the level density parameter a = 15.3 MeV~^ is also plotted with a dashed 
curve for comparison. 

Fig. 5 The calculated rotational E2 strength distribution as a function of the 7-ray energy 
E^. The four flgures correspond to the transitions decaying from the flrst, 6-th, 
24-th, and 200-th levels with (-, 1) in ^^^Yb at a; = 0.5 MeV, respectively. The 
magnitude of the strength is given as the fraction relative to the total strength. 

Fig. 6 The rotational strength distributions are shown for transitions from 50 levels av- 
eraged in each bin. The left-top quarter is for the flrst bin (from the flrst to the 
50-th level) with (— , 1). The other three quarters are for the levels from 51-th to 
100-th, 101-th to 150-th and 151-th to 200-th. The excitation energy from the 
yrast line is indicated for each bin. The smoothed distribution function is drawn 
with a solid curve. The transitions shown are for (— , 1) in ^^^Yb at a; = 0.5 MeV. 
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Fig. 7 The distribution of the normahzed rotational strengths for -"^^^Yb at a; = 0.5 
MeV. For convenience we draw the distribution as a histogram against the square- 
root yfsij of the normahzed strength. In this representation the Porter-Thomas 
distribution, shown as dashed curves, exhibits Gaussian shape. 

Fig. 8 The distribution of the normahzed rotational strengths. In this figure, the tran- 
sitions employed in the right-bottom quarter of Fig. 7 is divided into two groups 
according to whether their transition 7-energies are in the central part of the 
rotational strength function or in the tail part. 

Fig. 9 The branching number is shown as a function of the excitation energy from the 
yrast. The histogram is the result of the microscopic calculation while the dashed 
curve represents the Porter- Thomas limit given by Eq.(12) in which the rotational 
damping width Trot of 200 keV is adopted. The values are calculated for ^^^Yb at 
uj = 0.5 MeV. 

Fig. 10 The A3 statistics of the energy level fluctuation associated with the routhian 
eigenvalues e^(a;) in ^^^Yb at = 0.5 MeV. The eigenvalues are grouped into 
four bins in the same way as in Figs. 6 and 7. The long-dashed and the dashed 
curves represent the limits of the Poisson distribution and the GOE random matrix 
model, respectively. The fitted level density for (— , 1) is shown in the inset with 
the adjusted parameters of the constant temperature formula. 

Fig. 11 The E2 strength distributions as a function of the 7-ray energy for the lowest 9 
levels with (tt, a) = (+, 1) in ^^®Yb at u = 0.5 MeV. The branching number as 
well as the energy above yrast are also shown for each level. 

Fig. 12 The distribution of the normalized rotational strengths Sy, where the lowest 50 
levels for each (tt, a) are grouped into smaller bins of 10 levels. The approxi- 
mate energy regions corresponding to the bins are also shown. In this figure the 
transitions in the 10 nuclei around -"^^^Yb are sampled. 

Fig. 13 Distribution of the normalized nearest neighbor level spacings S for the unmixed 
rotational-band levels, defined by nbranch < 2, near the yrast at a; = 0.5 MeV (the 
upper half). The result for the independent-particle cranked shell model without 
the residual interaction at the same cu is shown in the lower half. The spacings are 
sampled from the ten nuclei around ^^^Yb and are normalized to give < 5" >= 1. 
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